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. . Abstract 

Test particles in a two dimensional, turbulent MHD simulation are found 
to undergo significant acceleration. The magnetic field configuration is a 
periodic sheet , pinch which undergoes reconnection. The test particles are 
trapped in the reconnection region for times of order an Alfven transit time 
in the large electric, fields, that characterize the turbulent : reconnection 
process at ...the relatively .. .large . magnetic .Reynolds number used in the 
simulation. The maximum speed attained by these particles is consistent with 
an analytic estimate which depends on the reconnection electric field, the 
Alfven speed,, and the ratio of Larmor period to the Alfven transit time. 
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Acceleration tOirfhigh ■uneTgies’uf a small fraction of the particles in a 
plasma is familiar^ in space physics, astrophysics and fusion physics. Two 

examples are the production of energetic particles in solar flares and in 
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planetary magnetospheres during magnetic substorms. Giovanelli , Dungey , 
Sweet 3 , Parker “ and others suggested that particles can be accelerated to 
high energies by strong electric fields which might arise near X-type 
u magnetic neutral points.. The question of whether strong electric fields 
'Sexist near such neutral points has led to the study of magnetohydrodynamic 
(MHD) reconnection. Conclusions about the intensity of the self-consistent 
•; electric field- near the X-poirrt and therefore, the ’’rate of reconnection" 
vary considerably in reconnection models’. The collisional tearing 
instability *^ Parker-Sweet mechanism 3,7 and Petschek's model 8 predict values 
of reconnection rates which decrease at high conductivity. 

The ^effects of turbulence in analytical and computational models have 
often been excluded by assuming that the background magnetic field is smooth 
initially and that perturbations are symmetric and/or infinitesimal. There 
is indication 9 ’ 10 ’ 11 that finite-am-plitude fluctuations can lead to 
turbulence in X-point dynamics , 'thereby elevating reconnection rates. 
Consequently, turbulent reconnection 'might maintain strong X-point electric 
fields (and the possibility of' significant particle acceleration) at high 
conductivity. 

The question of whether a subpopulation of plasma particles might be 
effectively accelerated by the reconnection zone electric field has also been 
extensively studied 13 17 . A typical approach has been to calculate the 
trajectories of "test-particles" in model MHD magnetic and electric fields. 
The electric field has usually been included parametrically, but there has 
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been at least one study of test particle orbits in the dynamic electric 
fields of a symmetric, forced non-turbulent MED simulation. It generally has 
been found that test particles do not spend very much time near X-type 
neutral points, but can easily be trapped near O-type neutral points 12 ’ 14 . 
This is a crucial point in the assessment of the importance of reconnection 
as a particle accelerator. In nonsteady incompressible simulations, strong 
accelerating**^ eleetr-re- fields do not appear near O-points which are the 
centers about which reconnected flux gathers after emerging from the 
reconnection zone. A long residence time near this type of ' O-point is 
unlikely to produce significant particle acceleration. The short residence 
times of^. particles- near the X-point region has been viewed as a limitation on 
the ef fic±encyt-^of-- ffre^X-point acceleration mechanism. 

Ih^tffTs letter, we describe an analysis of particle orbit tracing in MHD 
fields ~ which evolve and reconnect in the presence of finite amplitude 
fluctuations. We find that turbulent fluctuations appearing near the 
recottaection zone can trap test particles in the strong electric field region 
for long enough times to produce significant acceleration. Consequently the 
turbulent X-polnt mechanism is an efficient accelerator. 

The: fields in which the test particles propagate are produced by a 
spectral-method, incompressible MHD simulation technique that has been 
described in detail elsewhere 9 ’ 10 . The magnetic field B and fluid velocity 
field v lie in the x,y plane and vary in that plane, all quantities being 
independent of the third coordinate, z. Thus, the electric current density, 
magnetic vector potential, fluid vorticity and MHD electric field each have 
only single components in the z direction, designated J, a, <*> and E respec- 
tively. The familiar dimensionless units 19 * 9 include an arbitrary length 
scale L, a characteristic Alfven speed, V , and transit time of unit 

3 
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- distance, 'c —.X/V • 

' cl 

A.. 128*.;: spectral:. truncation and ; ant explicit ,: times tep of 1/512 were used. 
The trims .lasted 10&Alf^en;..£±mes» ..;iJ&e-..d±mensioni-ess. ..^fliiid and -magnetic 
Reynolds numbers (reciprocals of the dimensionless, viscosity, v and 
resistivity, W) were each equal to 1000 in these simulations. The values of 
u (k,t) and .a(k,t) for all values of 2— d wavevectors retained were stored 
every -At » 10/512 Alfven times. " The reconnection process is modeled as an 
unforced, turb ulen t initial value, problem, in which at t ** 0 magnetic . energy 
: resides , almost entirely in . a small, number of Fourier modes that describe two 
-ioppositely . directed neutral sheets. ... Turbulent ..reconnection.. is triggered by 
. adding af t = 0 a small but . finite -level of broad-band random-phased ."noise" 
perturbations. 

Test particles were ; introduced into the simulation fields at t = 0. 
Their position X and three dimensional velocity U were then advanced in time 
and space according to dX/dt “ U. and 

dU/dt « a(U x b + E) (!) 

where a = £it, fl = ZeB/mc r is.'fhe nominal test particle- gyrofrequency, and E = 
-v * B + id is the MHD electric field. 

* 

Eq. (1),. written with particle velocity in units of Alfven speed, 
emphasizes the role of the constant a in coupling the natural units of 
particle motion to those of MHD turbulence. Particle acceleration is 
entirely due to the z-component of electric field. The implementation of the 
particle-pushing algorithm uses the fact that the z-component of Eq. (1) 
implies 
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U z + a a = constant (2) 

This is an expression of the constancy of the z-component of canonical 
momentum and allows the effect of the electric field to be incorporated 
implicitly through the space and time dependence of the vector potential. To 
ensure that the particles sample the fields in the detail afforded by the 
high resolution MHD simulation, we employed a 2-d cubic spline interpolation 
procedure to calculate field values at particle positions. The fields were 
also linearly interpolated in .time, hence the electric field is constant 
between independent snapshots of the MHD simulation. Time integration of the 
x and y components of Eq. (1) was done by a fourth-order Runge-Kutta method. 
The times teps were chosen to be much less than the time between frames of the 
MHD fields to give good resolution of the particle gyromotion. 

The first run we describe consists of 500 particles randomly distributed 
in space, each moving with the unit Alfven speed at 45° pitch angle, but with 
randomly chosen gyrophase. The value of a is 643, so that the characteristic 
gyroradius is much less than the unit cell of the MHD simulation. Fig. 1 
shows magnetic field line plots in the area around the lower .of the two 
current sheets at several times. New magnetic islands are appearing due to 
the reconnection process. Furthermore, since the same field lines are 
plotted in each frame (same values of a(x,y)), it is clear that magnetic flux 
is being removed from the strong field regions above and below the current 
sheet . 

Figure 1 also shows the position of three of the test particles super- 
imposed on the field lines at each time. These particles were selected 
because of their significant energy gain, and because their trajectories were 
typical of particles gaining energy along either the upper or lower sheets. 
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Throughout the run, most of the particles are undergoing relatively 
- unperturbed gyromot ion ..in the strong magnetic afield regions away from the 
reconnection .zones . -nea r .vthe upper ^and . ^JLo ver^ -X— points^ . asJowe ver,. . Figure 2 
shows that the. particle energy distribution changes considerably during the 
run in that a small number of particles, including the three in Fig. 1, are 
accelerated to thousands of times their initial energy. It is apparent in 
Fig. ■ 1 that the particles were tra pped within small "magnetic bubbles” . in 
the reconnection .zone which are turbulent fluctuations of the. magnetic field. 

• Detailed : analysis . of ~the MHD simulation' 11 .has suggested that these 
... fluctuations- are -characteristic of reconnection triggered by finite amplitude 
=.. turbulence. • . The bubbles produce .multiple X— points within the reconnection 
aone which remains centered about a single filament of electric current 
density. Once formed, these bubbles move outward from the reconnection zone 
with the exhaust velocity of the fluid jets (“0.85V ). 

cL 

The . trapped particles move with the bubbles , cross the strong electric 
field region where they are accelerated, and enter the reconnected island 
region around an 0-point. Within the large island the particles experience 
less acceleration because the reconnection electric field is weak there. 
Some particles lose ^energy in the- 0— point region, presumably due to 
encounters . with field lines?, which-,. precede ; as “the ^reconnected island grows. 
This behavior.; is illustrated in Jig. 3 by the time history of the kinetic 
energy of the particles identified in Figure 1. 

The physics of the particle trapping and acceleration is easily under- 
stood. Each small magnetic." bubble is a trapping center for test particles. 
In this simulation, particles in the reconnection zone experience a fluctuat- 
ing electric field with average magnitude of £ “ 0.1 (in units of B V ). 
Particles trapped in a., bubble are entrained in the high-E region for about 
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one Alfven transit time. 

From the conservation law (Eq. 2) it follows that if a particle spends a 
unit time in a reconnection electric field, e, its velocity will change by 
AV z - ea. With <* = 643, the upper limit to the velocity is about 64. Note 
that in Fig. 3, particle no. 382 reaches a speed of about 60. This is the 
highest speed that we have found for particles in this simulation. A 
preliminary test of this scaling was performed using a second set of 500 
particles with a = 321.5. The highest speed attained was about 30. 

Assuming the validity of this scaling and neglecting the initial speed of 
the test particles, the maximum kinetic energy is e max 58 m(&V z ) 2 /2 = e 2 ct 2 m/2 
where e is in units of V^ 2 . This is equivalent to 

£ max " °' 5 mV a 2 t e a<- p /c)(Z/Z o )(m o /m)] 2 (3) 

where Z q and m Q are the plasma frequency, charge state and ion mass of 

the background plasma. Eq. (3) assumes that particle motion in the 
z-direction is unbounded. If the acceleration acts only until the particles 
have traveled a distance L in the z-direction, then the estimated duration of 
the acceleration process becomes t* = ^(2 /(eo)), and the maximum energy is 

‘.ax ' "a' t'tt« p /c)(z/Z 0 )(» 0 /*)] 

Note that this is proportional to the charge of the test particles, but 
independent of their mass. 

These simulations indicate that turbulent fluctuations in a reconnecting 
magnetofluid can trap test particles in the strong electric field region and 
accelerate them to speeds of order cSJtV . In a future publication, several 

Si 
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issues will be studied including the range of for which trapping occurs, 
the dependence of the magnitude - of 'the acceleration on Reynolds number and 
initial :.particle speeds,- and the . inf luence of a uniform DC magnetic field in 
the z-direction. 
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Figure Captions 

Figure 1. Magnetic field line plots in the area of the lower current sheet 
at several times. The positions of three particles are indicated by the 
symbols: the triangle is particle number 382, the circle is particle 358, 
and the square is particle 281. The contours label the same vector potential 
values in each panel. 

Figure 2. Particle energy distribution at t = 9, indicating' that most 
particles have gained little energy; but a small high energy population has 
been formed . 


Figure 3. The time history of the kinetic energy of particles 382, 358, 281 
(cf. Figure 1). 
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